|gr-qc/0505TT5| 
UFIFT-QG-05-02 
CRETE-05-11 



STOCHASTIC QUANTUM GRAVITATIONAL INFLATION 

in" 

O 1 N. C. Tsamis t 

O . 

CN . Department of Physics, University of Crete 

GR-710 03 Heraklion, HELLAS. 

R. P. Woodard* 

Department of Physics, University of Florida 
Gainesville, FL 32611, UNITED STATES. 



(N 
> 

in 



ABSTRACT 

in 

During inflation explicit perturbative computations of quantum field theo- 
ries which contain massless, non-conformal fields exhibit secular effects that 
grow as powers of the logarithm of the inflationary scale factor. Starobinskh's 
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q-i technique of stochastic inflation not only reproduces the leading infrared log- 

\ arithms at each order in perturbation theory, it can sometimes be summed 

to reveal what happens when inflation has proceeded so long that the large 
logarithms overwhelm even very small coupling constants. It is thus a cosmo- 
^ ■ logical analogue of what the renormalization group does for the ultraviolet 

logarithms of quantum field theory, and generalizing this technique to quan- 
tum gravity is a problem of great importance. There are two significant 
differences between gravity and the scalar models for which stochastic for- 
mulations have so far been given: derivative interactions and the presence of 
constrained fields. We use explicit perturbative computations in two simple 
scalar models to infer a set of rules for stochastically formulating theories 
with these features. 
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1 Introduction 



During inflation, the energy-time uncertainty principle allows any massless 
virtual particle that emerged from the vacuum to persist forever. If classical 
conformal invariance is present, the rate of emergence redshifts so that very 
few such virtual particles are present. It is only for gravitons and minimally 
coupled scalars that classical conformal invariance is broken in such a way 
that inflation can give strong enhancements of quantum effects. 

• de Sitter Inflation: A locally de Sitter geometry provides the simplest 
paradigm for inflation. To see why, consider a general homogeneous, isotropic 
and spatially flat geometry: 

ds 2 = — dt 2 + a 2 (t)dx-dx . (1) 

Derivatives of the scale factor a(t) give the Hubble parameter H(t) and the 
deceleration parameter q{t): 

-r-r, n o , . ad ^ H , . 

H(t) S - , q (t) s -_ = -!-_ . (2) 

The nonzero components of the Riemann tensor are: 

R° i0j = -qH 2 9ij , R M = H 2 (5\g 3i -5\g 3k ) . (3) 

Inflation is defined as positive expansion (H(t) > 0) with negative decelera- 
tion (q(t) < 0). On the other hand, stability - in the form of the weak energy 
condition - implies q(t) > —1. At the limit of q = — 1 we see from (JHJ) that 
the Riemann tensor assumes the locally de Sitter form: 



hm Rf alw = H 2 (5^ g uv - S" v g„») . (4) 

It follows from (0) that the Hubble parameter is actually constant, and that 
the zero of time can be chosen to make the scale factor take the simple 
exponential form we shall henceforth assume: 

de Sitter Inflation a(t) = e Ht . (5) 

• Particle Production: The homogeneity of spacetime expansion in ([Tfl 
does not change the fact that particles have constant wave vectors k, but it 
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does alter their physical meaning. In particular, the energy of a particle with 
mass m and wave number k becomes time dependent: 1 



E(t, k) 



k 2 



This results in an interesting change in the energy-time uncertainty principle 
which restricts how long a virtual pair of such particles with wave vectors 
±k can exist. If the pair was created at time t, it can last a time At given 
by the integral: 

rt+At 

dt' E(t', k) ~ 1 . (7) 

Just as in flat space, particles with the smallest masses persist longest. For 
the fully massless case, the integral is simple to evaluate, 



t+At 

dt' E(t', k) 





1 - e- HAt 


m=0 





(8) 



Ha(t) 



We, therefore, conclude that any massless virtual particle which happens to 
emerge from the vacuum can persist forever provided: 

Unbounded Lifetime ==> k < Ha(t) . (9) 

• Conformal Invariance: Most massless particles possess conformal in- 
variance. A simple change of variables defines a conformal time rj in terms of 
which the invariant element (0) is just a conformal factor times that of flat 
space: 

ds 2 = —dt 2 + a 2 (t)dx-dx = a 2 {rj) (—drj 2 + dx-dx^j , drj = —p- . (10) 

a(t) 

In the (77, x) coordinates, conformally invariant theories are locally identical 
to their flat space counterparts. The rate at which virtual particles emerge 
from the vacuum per unit conformal time must be the same constant - call 
it T - as in flat space. 2 

1 Of course "energy" is not a good quantum number in de Sitter background. What we 
mean by u E(t, k)" is the function whose integral times (— i) determines the phase of plane 
wave mode functions. This statement is exact for m — and true in the WKB limit for 
to ^ 0. 

2 For particles which do not possess conformal invariance, the rates are generically 
different in de Sitter than in flat spacetime ^ |2] ■ 
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Hence, the rate of emergence per unit physical time is: 

dN dN dr] T 
dt drj dt a(t) 



;n) 



Consequently - although any sufficiently long wavelength, massless and con- 
formally invariant particle emerging from the vacuum can persist forever 
during inflation - very few such particles will actually emerge. 

• Quantum Enhancement: Gravitons and minimally coupled scalars are 
two kinds of massless particles which do not possess conformal invariance. 
To see that - unlike massless conformally invariant particles - the production 
of these two kinds of particles is not suppressed during inflation, note that 
each polarization and wave number behaves like a harmonic oscillator with 
time dependent mass and frequency: 

11 k 
L = -mq 2 mwV , m(t) = a 3 (t) k uo{t) = —— . (12) 

2 2 Clitj 

The Heisenberg equation of motion can be exactly solved: 
k 2 

q + 3Hq + —q = =>- q(t) = u(t,k) a + u*(t,k) a t , (13) 
a 1 

where the mode functions u and the commutation relations obeyed by the 
operators a and oft are given by: 

H [ ik ( ik \ + 

u[t, k) = 1 — exp , [a , cr 



f 2k 3 



Ha{t) 



exp (— -| , [a.at] = 1 . (14) 



The co-moving energy operator for this system is: 

E(t) = \m(t)q 2 (t) + l -m{t)uj 2 {t)q 2 {t) . (15) 

Owing to the time dependent mass and frequency, there are no stationary 
states for this system. At any given time the minimum eigenstate of E(t) 
has energy \uj{t), but which state this is changes for each value of time. The 
state \ fl) which is annihilated by a has minimum energy in the distant past. 
The expectation value of the energy operator in its presence is, 



n 



E(t)\Q) = \a\t)\u{tM 2 + \a(t)k 2 \u(t,k)\ 2 = A + E^L . (16) 
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If we think of each particle as having energy A;a~ 1 (t), it follows that the 
number of particles N with any polarization and wave number k grows as 
the square of the inflationary scale factor: 



N(t,k) 



Ha{t) 
2k 



(17) 



Quantum field theoretic effects are driven by essentially classical physics 
operating in response to the source of virtual particles implied by quantiza- 
tion. On the basis of ()17j) . one might expect inflation to dramatically enhance 
quantum effects from massless, minimally coupled scalars and gravitons. This 
has been confirmed explicitly and the oldest results are the cosmological per- 
turbations induced by scalar inflatons [3] and by gravitons jlj. The more 
recent result which motivated the present analysis is that the gravitational 
back-reaction from the inflationary production of gravitons induces an ever 
greater slowing in the expansion rate [5J E] ■ 



• Quantum Cosmology: The Lagrangian is the two-parameter effective 
gravitational theory: 

C GR = -^—(-2A + R) ^9 , H 2 = \a>0. (18) 

What was actually computed [0] is the graviton one-point function, about a 
locally de Sitter background, in the presence of a state which is free Bunch- 
Davies vacuum at t = 0. However, if the resulting distortion of the back- 
ground geometry was viewed in terms of an effective energy density and 
pressure the perturbative infrared results would be: 

p(t) = sS? + 1 W^B" l2a + (H) + 0(k4) • (19) 
pW = -£g + (j ^ 1 {\^ a + ^} + °< k1) • < 20) 

where k 2 = lG^G is the loop counting parameter of quantum gravity. The 
one-loop effect from the kinematic energies of inflationary gravitons is a con- 
stant [3 |H] that must be subsumed into A. The next order effect is secular 
because the kinematic energies interact with the total graviton field strength 
which grows as more and more gravitons are produced. The induced energy 
density is negative because the gravitational interaction is attractive. 
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One would like to improve the aforementioned computation on the level 
of perturbation theory and beyond. 

• Invariant Regulation: A perturbative improvement would be to use 
.D- dimensional Feynman rules and dimensional regularization to obtain the 
fully renormalized answer and avoid complications that can arise from non- 
invariant counterterms. Such a procedure has been used in the simpler case 
of a massless, minimally coupled scalar field with a quartic self-interaction 
of strength A in a non-dynamical, locally de Sitter background: 

C v = -\yf=g gT d,tp - ^ V=g^ • (21) 

This theory can be renormalized so that, when released in free Bunch-Davies 
vacuum at t — 0, the energy density and pressure are P ITUj: 

+ 0(A 2 ) . 

(22) 
+ 0(A 2 ) . 
(23) 

This is a completely renormalized result that - besides the leading infrared 
term - explicitly exhibits all the sub-leading infrared pieces. The one-loop 
effect from the kinematic energies of inflationary scalars is again a constant 
that must be subsumed into A. The next order effect is secular because the 
ip 4 self-interaction involves the total scalar field strength which grows as more 
and more scalars are produced. The induced energy density is positive, for 
A > 0, because the tp A term adds to the energy density. In this model we 
might also guess that the effect is self-limiting because the classical restoring 
force tends to push the scalar back to towards zero energy density. 

• Infrared Logarithms: Both gravitation and the scalar model, as (jl9H20|) 
and (j22H23|) show, exhibit infrared logarithms - factors of m(a) = Ht. 3 As 
inflation proceeds, these infrared logarithms grow without bound until they 
eventually overcome the small coupling constants - (k,H) 2 for gravity and A 
for the scalar model. We cannot conclude that there is actually a significant 

3 Secular terms are ubiquitous in quantum field theory; for example see [iTHT2] . 
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Prcn(t) 



A 



XH 4 



8nG 2 6 vr 4 



1,2 2 . 
- In a H — a 
2 9 



is 

Z n=l 



n + 2 



-n-l 



n 



Pren{t) 



A 



XH 4 



3ttG 2 6 7T 4 



l l2 1 , 1 ^ n 2 - 4 
— In a H — In a H — > -, r- 

2 3 6-ejn + l) 2 



-n-l 



change in the background expansion rate because the higher order results 
remain unknown. The legitimate conclusions are rather that: 

(i) The expansion rate decreases for gravitation and increases for the scalar 
model; and 

(ii) Both effects eventually become non-perturbatively strong. 

In Section 2 we explain why infrared logarithms occur, both from the math- 
ematics of perturbation theory and on physical grounds. We also show that 
the field operator behaves like a stochastic random variable in the leading 
logarithm approximation. 

• Non-perturbative Extension: Since time evolution makes the secular 
growth of the infrared logarithms unbounded, perturbation theory eventu- 
ally breaks down. To reliably find out what happens after the breakdown, 
we must develop a non-perturbative technique. 4 This seemingly impossible 
task may have a reasonable chance to produce a satisfactory method be- 
cause inflationary evolution eventually makes the part of the quantum field 
responsible for particle production stochastic. Since it is precisely this part 
of the full field that drives the infrared effects of interest, we can try to iso- 
late its contribution in the equations of motion and solve for its evolution 
non-perturbatively. 

Starobinskii has long argued that his technique of stochastic inflation 
should recover the leading infrared logarithms at each order in perturbation 
theory. A known example where the leading infrared late time evolution 
has been calculated beyond perturbation theory using stochastic techniques, 
is the scalar model (|21j) : this was done by Starobinskii and Yokoyama [To] 
and its explication occupies Section 4. Their result agrees with the intuitive 
expectation that the growth of its field strength is eventually halted by the 
classical restoring force of the potential. 

The gravitational system (fTSj) to which we should like to apply the same 
method, differs in many ways from the particular scalar model (|21|). Two 
basic differences are the presence of derivative interactions and constrained 
fields in the effective gravity theory. Both differences must be addressed and 
the stochastic approximation rules must be extended so that they can be 
applied to (JTHJ). This is done in Section 3 where we also show that the field 
operator behaves like a stochastic random variable in the leading logarithm 
approximation. In Sections 5 and 6, we consider and analyze suitable scalar 

4 For non-perturbative methods different from the one we shall describe here see |13j . 
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models with derivative interactions and a reasonable analogue of constrained 
fields. Our conclusions comprise Section 7. 

The underlying idea behind the stochastic analysis of such theories is that 
the quantum field consists of two parts; a part which contains the ultraviolet 
effects that just redefine the parameters of the low energy interactions, and a 
part which contains the infrared effects from inflationary particle production 
and is responsible for the secular infrared logarithms. The ultraviolet sector 
decouples from the infrared except for: 

(i) The constant renormalizations of the low energy parameters its interacting 
component furnishes; and 

(ii) Its function as a reservoir of stochastic perturbations as the inflationary 
redshift pushes more and more modes into the infrared. 



2 The Physics of Infrared Logarithms 

The origin of the infrared logarithms in expressions (|19H2(J|) and (J22H23)) can 
be understood physically as well as from the mathematics of perturbation 
theory. Although the physical understanding is vastly more important in 
guiding the generalization we must make, we shall begin by explaining how 
perturbative computations are done invariantly in a locally de Sitter back- 
ground. Besides revealing the sources of the infrared logarithms, these tech- 
niques provide an invariant separation between infrared and ultraviolet de- 
grees of freedom. Moreover, these techniques shall be used in the computa- 
tions of Sections 5 and 6 which relate exact results from perturbation theory 
to stochastic realizations. 



• Invariant Regulation: We employ dimensional regularization in posi- 
tion space. The vertices are straightforward to obtain in any background, 
and the only difficulty comes in finding the propagators in arbitrary space- 
time dimension D. These propagators are expressed in terms of a de Sitter 
invariant length function we call y(x; x')\ 



y(x;x') = a(t)a(t') 



1 1 \- ! 

a a' 



(24) 



Its physical meaning in terms of the invariant length £(x; x') between x^ and 



x'^ is: 



y{x;x 



4 sin 2 



-H£(x;x') 



(25) 
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• Conformal Scalar Propagator: The simplest propagator is that of a 
massless, conformally coupled scalar [To] : 

*A CF (*;*') = ^r(^-l) (-f- 1 . (26) 
(471-)- ^2 > \y* 

Because y(x; x') vanishes when = x M , and because one always interprets 
D so that zero is raised to a positive power in dimensional regularization, 
the coincidence limit of the conformally coupled propagator vanishes: 

iA CF (x;x) = . (27) 

• Scalar Propagator: Although conformal invariance suppresses inter- 
esting quantum effects, the conformal scalar propagator is quite useful for 
expressing the propagator of the minimally coupled scalar EI] : 

iA A (x;x') = iA CF (x;x r ) (28) 
H D - 2 T(D-l) I D T 2 (f) ,4x#-a / 7t \ w ,.\ 

H D ~ 2 ~ J 1 T(n + D-1) f y^n_ 1 r(w + f + 1) f y^ n -%+2 



(47r)T n tU™ r (™ + f) V4y n-f + 2 I>+2) V 4 

This expression might seem daunting but it is actually simple to use in low 
order computations because the infinite sum on the final line vanishes in 
D = 4, and each term in the series goes like a positive power of y(x;x'). 
This means that the infinite sum can only contribute when multiplied by a 
divergence, and even then only a small number of terms can contribute. 

• Correlation Source: The explicit factor of ln(aa') present in (J28|) is 
one source of infrared logarithms. It gives the secular dependence of the 
coincidence limit [TTHTSl IP9] : 



iA A (x;x) = N g ; \ -7rcot(-£>) + 2\na\ . (29) 



H D-2 Y(D-l) f /7T 
n 7JT\ — S COt — . 

(4tt)# T(f) I V 2 

This factor of ln(aa') is also interesting in that it breaks de Sitter invariance. 
Of course tp(x) transforms like a scalar; the breaking derives from the state in 
which the expectation value of the two free fields is taken. Allen and Folacci 
long ago proved that the massless, minimally coupled scalar fails to possess 
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normalizable de Sitter invariant states (20] • A final point is that the factor 
of ln(aa') is actually an invariant, it just depends upon the initial value sur- 
face at which the state is released. In fact, it is the Hubble constant times 
the sum of the invariant times from x M and x' M to this initial value surface. 
The physical reason for the appearance of such a term is the increasing field 
amplitude due to inflationary particle production. 

• Graviton Propagator: We define the graviton field ijjp^( follows: 

9tiu(x) = a 2 (77^ + ' K 2 = lQ-n:G . (30) 

Our gauge fixing Lagrangian takes the form: 



-^a D - 2 V ^F„F, 



p ± V 1 



(31) 



where F M is an analogue of the de Donder gauge fixing term of flat space [2T] : 



pp,a 



^/W + (D-2)Ha^ p 5° a ) . 



(32) 



Because space and time components are treated differently it is useful to have 
an expression for the purely spatial part of the Minkowski metric: 



V = ^ + <? 6 u 



(33) 



With these definitions the graviton propagator takes the form of a sum of 
three constant index factors times three scalar propagators: 



pu^pa 



(34) 



I=A,B,C 

The explicit expressions for the index factors are |2*2*] : 

2 



aw p<j 

rpB 

liv L per 

rpC 

li-v 1 pa 



2 V p,(p V 'a)u JJ_r^VpuVp<T ) 



'{p'lv){p u a) 5 

2 



(D-2)(D-3) 



(D-3)5X + V + 



I pa 



(35) 
(36) 
(37) 
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The A-type propagator is identical to J2Ej), while the B-type and C-type are 
given by: 



iA 



B (x;x 



iA r 



[x: x 




(38) 



iA c (x\x') = iA CF (x;x') + _ 

(4?rJ 2 n=0 l 

/ d n r 

-In 1-31 - 

T(n + 2) 

For completeness we also give the ghost propagator in this gauge: 



,A, 



(x;x') = r]^ iA A (x;x') 



- 5j<£iA fl (x;x' 



(39) 



(40) 



Note that the infinite sums in ("38|) and ([3*9"]) vanish for D = 4 so that, in 
this limit, the S-type and C type propagators both agree with the confor- 
mal propagator. It is significant that only the A-type propagator contributes 
infrared logarithms. Since only physical graviton modes can experience the 
inflationary particle production responsible for infrared logarithms, we ex- 
pect only the A-type propagator to contain them and this is indeed the case. 
A final point is that the breaking of de Sitter invariance apparent in these 
infrared logarithms is a real effect, not an artifact of having employed a de 
Sitter non-invariant gauge. One way to prove this is by defining a de Sitter 
transformation of the graviton field to include the compensating diffeomor- 
phism needed to restore the gauge condition *2*3*j . 

• Volume Source: Although any infrared logarithms which appear in 
one-loop diagrams can only have come from the A-type propagator, there 
is another mechanism that can produce infrared logarithms in higher loop 
results such as (|19H20|) and (j22H23|) . This other mechanism is the growth of 
the invariant volume of the past light-cone from the observation point back 
to the initial value surface: 



dt' a' 



2^ 



H 



-D 



(£>-i)r. 



D-l- 



y(x;x' 
In a + 0(1) 



(41) 
(42) 
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Factors of this quantity arise naturally whenever undifferentiated propagators 
connect an interaction vertex - at x' M - with the expectation value of some 
observable - at x^. 

To obtain true expectation values for cases - such as cosmology - in 
which the "in" and "out" vacua either do not agree or are not even well 
defined [21], the Schwinger-Keldysh formalism [2H] must be employed. In 
this formalism the only net corrections come from interaction vertices which 
lie on or within the past light-cone of some observation point. The more 
familiar "in- out" matrix elements would harbor virulent infrared divergences 
from integrating over the exponentially large inflationary future volume [26 . 
The causality of the Schwinger-Keldysh formalism regulates these infrared 
divergences but simple correspondence implies that the regulated expressions 
must grow without bound at late times. 

3 Infrared Dynamics and Their Rules 

A crucially important consequence of our use of an invariant ultraviolet reg- 
ularization is that there should be a constant dynamical impact from all 
modes whose physical wavelength ranges from zero - the far ultraviolet - to 
any fixed value. Recall that quanta are labeled by constant wave vectors k, 
and that the mode with wavenumber k = \\k\\ begins to experience significant 
inflationary particle production when the number of particles N(t, k) > 1 or, 
equivalently, when k < Ha(t). This suggests the following physical separa- 
tion between "infrared" and "ultraviolet" modes: 

Infrared H <k< Ha{t) , (43) 

Ultraviolet k > Ha{t) . (44) 

The range of the physical wavelengths A p h = 2nk~ l a(t) of the ultraviolet 
modes is from zero to the constant if -1 . With any invariant regularization, 
the dynamical impact of such modes must be constant because they lie in an 
invariantly defined range. 

Had we taken the lower limit in expression (|43|) down to k — 0, the 
infrared phase space would also have extended over a constant physical range. 
However, taking k down to zero has long been known to result in infrared 
divergences [23 • We regulate these by working on the manifold T D ~ l x 3?, 
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with the range of the toroidal coordinates equal to a Hubble length |2o] : 

- <x i < ^H- 1 . (45) 



Other regulating techniques exist [28J but they all cause the effective infrared 
phase space to increase as the universe inflates. Indeed, this growth is the 
physical source of infrared logarithms. What happens is that the average 
field strength increases as it receives contributions from more and more in- 
frared modes. If any interactions involve the undifferentiated field, then this 
growth can enter into physical quantities. 

• Self-interacting Scalar: These assertions can be confirmed in the context 
of a self-interacting scalar theory in D = 3 + 1 dimensions and in the presence 
of the inflationary background (J5J): 

c = -\^g gT d,# d v <p - V{<p) , (46) 

where, for stability reasons, the potential V(<p) is to be bounded from below. 
The resulting field equation is: 

(p + 3H (p - — VV + V%) = . (47) 



The full - ultraviolet plus infrared - perturbative initial value solution can 
be obtained by iterating a Yang-Feldman equation [2H] for which the 
time has been set to t = 0: 



"m" 



<p(t,x) = cpa(t,x)-j dt' a 3 (t') j d 3 x' G Tet (t, x ;t' ', f) V'(tp)(t' , x') . (48) 
The retarded Green's function equals: 



(49) 



v ' 47r [ aa'H\\x — x'\\ 

+ 9(h\\x-x'\\ + --^)) . 
v a a ' J 

The free field ip (t,x) is expanded in terms of mode functions u(t,k) and 
operators a(k) and a\k) obeying canonical commutation relations: 

<p (t,x) = [ \ e its u{t,k) a{k) + e~ it3 u*{t,k) a ] (k)\ , (50) 
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where: 



u(t, k) 



H 



1 - 



ik 



ik 



a{k), a\k')] = (2n) 3 5 3 (k-k') . 



(51) 
(52) 



It should be noted that although (po(t, x) is only the lowest order part of the 
solution, it and its first derivative agree exactly with the full field on the 
initial value surface. 



• Infrared Field: To excise the ultraviolet modes (J44j) . we iterate what is 
essentially the same equation: 

$(t,f) = $ (t,f) - / dt' a 3 (t') fd 3 x' G iet (t,x ;t',x') V '(<$>) (t' ,x) , (53) 



but with the zeroth order solution restricted to only infrared modes: 



$o(t, x) = J -0^ 9(Ha(t) - k) je <&2 «(t, k) a % + e - ilg u*(t, k) atj . (54) 



This model is completely free of ultraviolet divergences, so we are justified 
in taking D = 3 + 1. To see that the model also reproduces the leading 
infrared logarithms at tree order it suffices to take the vacuum expectation 



value of 5 



$l(t,x)) = ( J^-e(Ha-k)^{l 
0V ' ' /vev J (2tt) 3 V J 2k 3 \ 



H 2 a 2 



H 2 f Ha dk { k 2 
1 + 



Ait 2 Jh k 

2 ( 



H_ 
2^ 



In a 



H 2 a? 



2a 2 



(55) 
(56) 
(57) 



Comparison with (j2H|) for D = 3 + 1 reveals exact agreement between the 
ln(a) terms. 



5 Vacuum expectation values are taken throughout in the presence of the free Bunch- 
Davies vacuum at t = 0. 
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• Infrared Field Equation: Our purely infrared field operator $(£, x) does 
not quite obey the original field equation (jlT|) because the kinetic operator 
fails to annihilate <&o(t,x). The reason is that, upon taking time derivatives 
of &o(t,x), there are extra contributions due to the presence of the time 
dependent upper limit Ha(t) of the mode sum in ([54)1 . Thus, the action of 
the kinetic operator on <& (t,x) produces momentum space surface terms: 

... V 2 

$ (t,x^ + 3H<S> (t,x)-—-r<$>o(t,x) = J r (t,x) + G(t,x)+3HJ r (t,x) , (58) 



a 2 (t) 



involving the sources: 

_ . ., , i (1 k' 

F(t, x) 

G(t,x) 



Hf ——k5(k-Ha{t))\e ils u{t,k)a % + e~^V(t, k) a\\ (59) 



H 



(2tt) 3 
d?k 
(2^F 



kS[k-Ha{t))le ik - s ii{t,k)a R + e~ lk - £ u*(t, k) aU (60) 



ik-x 



Note that the mode functions and their derivatives are simply constants at 
k = Ha(t): 

H 2 



u(t, k) 
u(t, k) 



H 



k=Ha{t) 



H 2 



l-i)e 



\u(t,k)\ 2 
\u(t,k)\ 2 



k=Ha{t) 



k=Ha(t) 



k 3 

El 



(61) 
(62) 



k=Ha(t) 

In view of (J5SJ) the equation obeyed by the infrared field <£>(£, x) is not 
(HZJ) but rather: 



3if $ 



o . 



(63) 



• Infrared Conservation: The infrared field equation ()63)1 does not leave 
the original stress-energy tensor conserved because stress-energy is being con- 
tinually dumped into the truncated system by ultraviolet modes which red- 
shift past the horizon. We can account for this by modifying what we call 
T^ v . To motivate the modification it is useful to write down the form we 
expect for the divergence of the stress-energy. For the zero component we 
should get: 



rp \H 

1 0fj, 



—Too — H 



3T o + <f%' +9 ij T 0i ,, 



(<£ - T - Q) + 3H($ -J 7 )- — $ + V'{$) 



(64) 
(65) 
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while for the spatial components: 



Toi — 3HToi H — -Tj 



(66) 
(67) 



We can enforce (J65|) and (|67|) with a stress-energy of the form: 
T 00 = ^($-^) 2 + -^V$-V$ + V(<l>) , 



T, 



Qi 



$-^) 2 + ^V$-V$ + V($) 



5 



(68) 
(69) 

(70) 



Once this form is assumed, the non-local source terms of the purely spatial 
components are determined by conservation: 



St 
S L 

s 



a 

a 2 d k 
V 2 

2 



Gd& + {®-T)dir 



Qd h <5> + {®-F)d k F 



a 4 did k 



a 



$- < F-0 + 3i7($-.F)- 



gd k $ + {$-F)d k T , (71) 

(72) 

{§-F)Q -HVF-V® (73) 

(74) 



a 

77 



We have, therefore, achieved a completely consistent model of just the in- 
frared modes that reproduces the leading infrared logarithms. Note that this 
establishes the decoupling of the ultraviolet sector. 

• Leading Infrared Field: Although the zeroth order field $o(^> %) contains 
only infrared modes, it is still quantum mechanical in that the field and its 
first time derivative do not commute: 



d 3 k 



6{Ha{t)-k 



uu — u u 



ik-Ax 



(75) 
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4ir 2 a 3 



11 " 2 sin (fc Arc) 



dk k l 



H 



A7r 2 (aAx) 



kAx 

sm(aHAx) — aHAx cos(aHAx) 
- sin(if Ax) + H Ax cos(i/Ax) 



(76) 
(77) 



where Ax = x — x! and Ax = ||Ax||. However, the leading infrared logarithm 
in (J57)) derives entirely from the constant first term of the long wavelength 
expansion of the mode function: 



u(t, k) 



H 



/2k 3 
H 

7W 



ik 
1/ k 



i k 



1 ( 



A: 



+ 0(k 4 



(78) 
(79) 



By constructing a free field, which we shall call <fto(t,x), with only this first 
term as its mode function: 



Leading Infrared 



k) 



H 



'2k 3 



(80) 



we would get precisely the same leading infrared logarithm: 



4>o(t,x) = 



d 3 k „/ . , , \ H 
(2^ ( Ha{t) - " 



Ak-x 



2k 3 



olz + e 



[81) 



Because - unlike a usual quantum field - the creation and annihilation parts 
of (po(t,x) are both multiplied by a phase factor with identical time depen- 
dence, the leading infrared field commutes with its time derivative: 



[(f>o{t,x) , (f) {t',x') 



. 



(82) 



Consequently, 0o behaves like a classical variable except for the operators 
a and which can take random values. Such a random but commuting 
variable might be termed "stochastic" . 



• Leading Infrared Field Equation: Replacing $o with 0o i n the basic 
equation (jo^j) does capture the leading infrared logarithms, but the resulting 
field theory is still more complicated than necessary. The leading infrared 
logarithms are preserved if we retain only the first term in the long wavelength 



16 



expansion of the retarded Green's function. From the expression: 



G ret ( t, x ; t', x f ) = iO(t- f) x (83) 



d 3 k 



u(t, k) u*(t', k) - u*(t, k) u(t', k) 



ik-Ax 



we can determine the long wavelength expansion of the bracketed term using 
dZOJ): 

u(t, k) u*(t', k) - u*(t, k) u(t', fc) ] = - 1) + Q(k 2 ) . (84) 



Retaining only the leading term from (jHljl in (jHojl we conclude: 

G iet (t,x ;t',x > ) — -Le( t -t')5 3 (x-x<) (^-^J • (85) 

Recall that the Green's function is multiplied by a factor of a' 3 from the 
measure of integration in the Yang-Feldman equation. Whereas the term 
(a' /a') 3 = 1 contributes coherently over the full range of integration, the 
term — (a'/a) 3 is only significant for t' ~ t. Hence, the second term in (I85J) 
is irrelevant in the leading logarithm approximation: 

Leading Infrared =^ a' 3 G™ t (t, x ; t', x' ) ~ — 9{t—t') 5 3 (x-x' ) . 

(86) 

We are at length lead to the following simple iterative equation for recovering 
the leading infrared logarithms: 

<f>(t,x) = (j> Q {t,x) - ^Lj*dt'V'(<f>(t',x)) . (87) 

To infer the local differential equation which the full stochastic field (f> 
obeys, we note that the time derivative of the stochastic free field 0o is a 
momentum space surface term: 

r d 3 k H 2 f 1 

Mt,x) = / — - 8(k-Ha(t)) -= e ik -*a % + a\ , (88) 



(2tt) 3 V v V V2fc 

U{t,x) . (89) 
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Taking the time derivative of (|87|) gives the Langevin equation for the stochas- 
tic field 4>{t,x) first obtained by Starobinskii [T41 I15j: 

<j>(t,x) = U(t,x) - -^V'(4>{t,x)) . (90) 

Physically, f^ is a stochastic source caused by the ultraviolet modes that 
are instantaneously redshifting across the horizon at time t. Because each 
increment is uncorrelated, this source represents Gaussian white noise: 

( f v (t, x) U(t', x) ) vev = ^S(t- t') . (91) 

The sort of system comprised by (j89H90j) has been much studied [3U|. Expec- 
tation values of functionals of the stochastic field can be computed in terms 
of a probability density g(t, (f>) as follows: 

( F[(f)(t,x)} ) = / du g(t,u) F{u) . (92) 

\ / VEV J —OO 

The probability density satisfies a Fokker-Planck equation whose first term 
is given by the interaction on the right hand side of equation (j9*Uj) and whose 
second term is fixed by the normalization of the white noise (|91|) : 



1 d 



3H 



H 3 d 2 



= T777^t V'(cl>)g(t,4>) + QM) ■ ( 93 ) 



8tt 2 



• Leading Infrared Rules: Relation ()90|) derives from the original field 
equation (jTfj) by applying the following rules: 

(I) At each order in the field — ip , (p 2 , and so forth — retain 

only the term with the smallest number of derivatives. (94) 

(II) For the linear terms in the field, each time derivative has 

a stochastic source subtracted. (95) 

For the kinetic part of (jlZJ), by applying rule I we select the Hubble 
friction term and, thereafter, by applying rule II (|93j) a stochastic source 
must be subtracted off because we are at linear order: 6 

V 2 

(p + ?,H^-—p — > 3Hip — > ZH^-U). (96) 



6 Since the scale factor varies much faster than the field during inflation, it is preferable 
to have derivatives act on the scale factor instead of the field. Thus, the single time 
derivative of the Hubble friction term dominates over the second time and space derivatives 
terms. 
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There are no derivatives or linear terms in the interaction part of (|47|) so we 
merely replace the full field by its stochastic counterpart: 

V\ip) — V{4>) ■ (97) 
Hence the full equation (jTTj) reduces as follows, 

f + SHtp-^cp + V'fr) — > 3HU-U)+V'(c/>) , (98) 

which is indeed equivalent to (|9~U|). 

It is important to note that we would get the same leading infrared loga- 
rithms by solving ()63|) for $. The enormous advantage of solving the stochas- 
tic equation (}90|) instead, is that the field (ft can be regarded as classical. 
Whereas there is little hope of being able to exactly solve the Heisenberg 
field equations for an interacting, four dimensional quantum field, classical 
field equations can sometimes be solved. 



4 Scalar Field with Quartic Self-Interaction 



For the quartic self-interaction (}2~T]) the full field equation is: 

V 2 A 



(f + 3Hif -if + - <p = 

a z o 

Applying the two stochastic reduction rules gives: 

A ,„ 



. 



+ 



6 







(99) 



(100) 



• Non-perturbative Solution: To compute expectation values we need the 
probability density g(t, 0) which obeys the Fokker-Planck equation: 



d_ 

dt 



g(t,4>) 



A d 



18H 



3 g(t, 0) 



+ 



H 3 d 2 



In 2 d4> 2 



Q{t,4>) 



(101) 



To recover the non-perturbative late time solution of Starobinskii and Yoko- 
yama we make the ansatz: 



lim g(t,(p) = goo{4>) 

t — >oo 



(102) 



19 



because the scalar force should eventually balance the tendency of inflation- 
ary particle production to force the scalar up its potential. This ansatz results 
in a first order equation: 



£?oc(0) ~ 
The solution is straightforward: 



(103) 



4 



7T 2 / x ; 



(104) 



and its non-perturbative nature is clear. It follows that the stochastic ex- 
pectation value of the 2n-th power of the field has the following late time 
limit: 

lim ( 2n (t,f) ) = / du goo (uj) u 2n , (105) 

t^OO \ I VEV ./_™ 



r(f + |) (m 



4\ 2 



r(|) W 2 \ 



(106) 



• Perturbative Agreement: To make contact between stochastic expecta- 
tion values and perturbation theory, we first differentiate the vacuum expec- 
tation value of <p 2n and then we use the Fokker-Planck equation [3*T| I32j: 

— ( <p 2n (t,x) ) = duo g(t,u) u? n , (107) 



Qt \ 1 ' /VEV 

n(2n-l)H 3 , - _ 2 . . v nX 



4tt 2 \ r v ' ' vev g# 

This relation can be more revealingly expressed in terms of a new time vari- 
able a and a rescaled coupling constant A: 



2 



1 - 4tt 
a = — In a , A = A . (109) 

47T" 1 9 



A differential recursion relation emerges: 

/ / 6 \2n\ , / , <t> \2n-2\ — / / 6 -s2n+2 



VEV 



(110) 
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whose solution has the form: 

(j) \2n* 

H 



(2n-l)!! a n F{Xa 2 ) , 



111^ 



where the function F n (z) obeys: 

2zF' n {z) + nF n (z) = nF n ^(z) - n(2n + l)zF n+1 (z) . (112) 
It is straightforward to obtain the first few terms of the series expansion: 

F n (z) = l--(n + l)z + —(35n 3 + 170n 2 + 225n + 74)z 2 + O(z 3 ) . (113) 
2 2oU ^ ' 

Hence, the stochastic expectation value of <p 2n has the following time depen- 
dence in perturbation theory: 



(2n- 1)!! 



H 2 
4tt 



In a 



n , A , n 



+ -^-(35n 3 + I70n 2 + 225n + 74 
280 ^ 



A 2 12 

In a 



36tt 2 



+... (114) 



Part of this result is in precise agreement with detailed and highly non-trivial 
calculations of the same quantities using quantum field theory 



5 Derivative Interactions 

A basic difference between the scalar model of Section 4 and gravitation, is 
the presence of derivative interactions. Since our main purpose is to derive 
the proper stochastic equation for quantum gravity, we first study a simple 
scalar model with derivative interactions defined by the following Lagrangian: 



~ V=9 <T d,A d v A- l -^~g <T d,B d v B 
A 



-g A 2 d.B d v B 



(115) 



Whereas A and A have the same equal-time commutation relation as a free 
field, B and B do not: 



A(t,x) , A(t,x") 
B(t,x) , B(t,x') 



iS 3 (x — x') 
a 3 (t) ' 
iS 3 (x — x 1 ) 
a 3 (t)[l + \A*{t,x) 



(116) 
(117) 
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This means we will need to include a homogeneous term in writing the Yang- 
Feldman equation for B. 

The full equations of motion are: 

' |f = - ^gTAB^Br = , (118) 



/=g 5 A » 2 

' || = B^ + Xg^AA,,B >v + ~A 2 i^ = . (119) 



Each equation of motion can be re-written to give the d'Alembertian of the 
field in terms of lower derivatives: 

A^ = \AB„B* , (120) 

\AA U B>» . , 

B*„ = £ . 121 

1 + ^A 2 



The full perturbative initial value solution comes from iterating the Yang- 
Feldman equations: 7 

\ ft r 

A(t,x) = A (t,x) - - J dt' a 3 {t') J d 3 x' G ret (t,x ; t', f ) x (122) 

A{x')B^{x')B^{x') , 
B(t,x) = B {t,x) - -J d 3 x' G TCt (t,x ;0,x) A 2 {0,x) B{0,x) (123) 



+ A J dt' a 3 (t') J d 3 x' G TCt (t, x ; t', x' 



, ^ ( A(x') AJx') B^(x') 



l + ±A 2 {x> 



Because both A and B are massless, minimally coupled scalar fields, the 
retarded Green's functions in (jl22H123|) are identical to (}4*9*j) . The free field 
expansions are similarly identical to (j50H52|) . except that each field has an 
independent canonically normalized set of creation and annihilation opera- 
tors: 



A (t,x) = I ^3 \ e**u(t,k) a % + e**u\t,k) at \ , (124) 



B (t, x) = J** I e**u{t, k) /% + e-^u*(t, k) 0j, \ . (125) 



(2k) 



7 Field arguments are sometimes compressed to 4- vector notation: (t,x) = x. 
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The second homogeneous term in (|123|) accounts for the non-canonical 
commutation relation (|117j) . This second term can be expressed as the inte- 
gral of a total derivative: 

A 



- J d 3 x' G Tet (t, x ; 0, x* ) A 2 (0, x') B(0, x') 



A r f 



2 Jo 
A ft 



2 Jo 



dt' j d 3 x'd^ I yj-g(x') g^(x') G iet (x; x') A 2 {x') d' u B(x') 
dt! a! 3 ( d 3 x' d'^G rct (x]x') g^{x') A 2 {x') & v B{x') 



A j dt' a' 3 j d 3 x' G TCt (x; x') 



A(x') A^{x')B^(x') 
. l + ±A\x>) 



(126) 



In reaching this final expression we have made use of the equation of motion 
(I121j) . Combining terms gives our final form for the Yang-Feldman equation 
of the field B, 



B(t, x) 



B (t,x 

A rt 



(127) 



- ^ y dt' a' 3 J d 3 x' d'^et (t, x ; t, /) g^(x') A 2 {x') d' v B{x') 



• Stochastic Realization: We need to make the transition from the full 
fields A and B to the stochastic fields A m and B m . Most of this proceeds as 
the reduction of Section 3, with slight modifications to accomodate derivative 
interactions. The free fields suffer the same truncation of their ultraviolet 
modes but we must now include the next order term in the long wavelength 
expansion of the mode functions because the first term has zero time deriva- 
tive: 



Ao(t,x) — > A m0 (t,x) 

d 3 k / , N x /-/ 
9{Ha(t) - k 



(128) 



f2vr) 3 V 



2k 3 



1 

1 + - 



k 



2VF a (t) 



B {t,x) 



B m0 (t,x) 



d 3 k „/„ ts jN H 
9[Ha{t) - k 



(2vr 



2k 3 



k 



2^Ha(t) 



) 



ik-x 



Ak-x 



at + e 



-ik-s + 



(129) 



[3 + e -**$ 



The stochastic limit of the retarded Green's function is unchanged: 



a' 3 G Te ±[t,x :t'x 



ret \ l/ i Jj ) v , ■ 



— 9{t-t') 5 3 (x-x') 



(130) 
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The new feature of these equations is that each interaction contains two 
derivatives. From the analysis of Section 3 it is apparent that differentiating 
a propagator precludes it from contributing an infrared logarithm. Said 
another way, it is always preferable to have a derivative act on the rapidly 
varying scale factor rather than the slowly varying field. We, therefore, expect 
that the Hubble friction term dominates the scalar d'Alembertian: 

... V 2 

B^ u = -B - 3HB + — B — ► -3HB m . (131) 

There are no d'Alembertians in our Yang-Feldman equations but the time 
derivative of the retarded Green's function in (|127jl contributes a similar 
term. The space derivatives always give two derivatives acting upon fields, 
and no term is comparably important to (|131|) when the two derivatives act 
on different fields: 

B U B* = -B 2 + ^VB-VB — ► 0, (132) 

cr 

An B'** = -AB + ^VA-VB — ► 0. (133) 

er 

By dropping all double derivatives of the fields we arrive at the equations 
which recover the stochastic fields: 

A m (t,x) = A m0 (t,x) , (134) 

B m (t,x) = B m0 (t,x) - £ f dt' A 2 1R (t',x)B m (t',x) . (135) 

2 jo 

As in Section 3, the local stochastic field equations are obtained by dif- 
ferentiating the stochastic Yang-Feldman equations (jl34H135|) with respect 
to time. Also, as in Section 3, derivatives of free fields produce stochastic 
source terms: 

A 1K0 {t,x) = f A (t,x) (136) 



(2vr) 3 V v, ^ 

1 r d 3 k 



-ik-x t 

k 



a k + e ik '* a l 



a 2 (t) J (2tt) 3 

B m0 (t,x) = f B (t,x) (137) 
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= / Wf S(k - Ha{t)) Wk { **** + e ~"^ } 

Because no derivatives of Am appear in the Yang-Feldman stochastic field 
equations we can neglect the final term in (J136|) but it must be retained 
in (|137|) on account of the B m which appears in (|135|) . The resulting local 
stochastic equations are: 

A m {t,x) = f A (t,x) , (138) 
B m (t,x) = f B (t,x) - ^Al(t,x) B m {t,x) . (139) 

These are in good agreement with the equations which are obtained by di- 
rectly applying the reduction rules I-II (|94H95|) of Section 3 to the full equa- 
tions of motion (111811 191) : 

-3H(A m -f A ) = , (140) 
-3H(B m -f B ) - ^A 2 m B m = . (141) 



Note the curious fact that interactions containing derivatives are free of 
stochastic source terms. 

It is amusing to note that we can obtain explicit operator solutions to 
the stochastic field equations (jl38H139|) . The exact solution for A m follows 
trivially from ()138j) : 

A m (t,x) = f dt' f A (t',x) = A m0 (t,x) , (142) 
Jo 

and shows that A m receives no corrections to its free field value A m0 . A few 
simple re-arrangements in (|139|) give the closed form operator solution for 

-Bib.: 



B ^ s) = L dt ' TTf^^) • (143) 

From ()143|) we conclude that B 1R exhibits a kind of spacetime-dependent 
field strength renormalization whereby each set of /5-modes which experi- 
ences horizon crossing is attenuated by the factor 1 + ^Af n0 . In other words, 
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Figure 1: The Feynman rules and the diagrammatic expansion to 0(A) of the 

vacuum expectation value of the squares of the fields in the scalar model 
with derivative interactions. Solid lines correspond to the scalar field A and segmented 
lines to the scalar field B. 
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the stochastic field B lK acquires corrections from its free field form which 
diminish with time and it reaches some constant value asymptotically. 



• Perturbative Agreement: It is possible to calculate the one and two loop 
contributions to the expectation values of the squares of the full fields A and 
B. Because they are both massless and minimally coupled fields, they have 
the same propagator given by (|28j) . The basic vertex is: 

V u ,2s(x) = -i^V^gg^dxvdto , (144) 

where we have used the notation of Figure 1 which also depicts all relevant 
graphs. In the Appendix we use these Feynman rules to compute: 

A 2 (t,x)) = (— )\na + 0(A 2 ) , (145) 

B 2 (t,x)) m = (^'ina + ^[-ln 2 « + O(lna)] + 0(A 2 ) . (146) 

It is apparent that ()145j) agrees with the stochastic result ()142j) . To see 
that ()146|) agrees as well with ()14Hj) . we perturbatively expand the stochastic 
operator solution: 

B m (t,x) = B m0 (t,x) - ~ f dt' A 2 m [t',x) / B (t',£) + 0(A 2 ) , (147) 

2 Jo 

and square it: 

B^(t,x) = B* n0 (t,x) - A B m (t, x) f dt' A 2 m {t\x) f B (t',x) + 0(A 2 ) . 

Jo 

(148) 

Taking the vacuum expectation value of ()148|) gives: 



5 2 R (t,x)) vEv = (E 2 R0 M)) VEV (149) 
-A f^dt' ( A 2 1R0 (t',x) ) yEv ( B m0 (t,x) Ut\x) ) vev + 0(A 2 ) . 
For the case in which t > t' we have: 



H 2 r dk ( 

B m0 (t,x)f B (t',x))_ = —j—\H 2 a'5{Ha'-k) (150) 

-— H Ha- -*)} = S5 . 
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Note that we discard the time dependent corrections to the wave function in 
B m which would in any case give sub-dominant terms like a~ 2 a' 2 . Hence we 
obtain: 



which establishes the non-trivial agreement between (|143J) and ([146)1 . 



Besides dynamical degrees of freedom, gravitation contains constrained fields. 
Because these fields possess no dynamical degrees of freedom there should be 
no independent stochastic sources for them. However, a constrained variable 
can depend non-linearly and non-locally upon the dynamical fields and it is 
by no means clear how to include the extra stochastic jitter the constrained 
field acquires due to each increment of dynamical modes which experience 
horizon crossing. 

No purely scalar model can mimic the gravitational system we really need 
to understand. It would be interesting - and it seems feasible - to gain insight 
into constrained degrees of freedom by attempting a stochastic formulation 
of scalar quantum electrodynamics. This model certainly possesses a gauge 
constraint, and the one loop vacuum polarization has been shown to contain 
an infrared logarithm |34j . However, as a first step we shall attempt to gain 
insight from purely scalar models by exploiting the existence of a covariant 
gauge for gravity in which all components of the perturbed metric are either 
minimally coupled or else conformally coupled 35 . 

The model we shall study consists of a massless, minimally coupled scalar 
A, a massless, conformally coupled scalar C and an elementary interaction 
between the two: 



VEV 




+ 0(A 2 ) 



(151) 



(152) 



6 Constrained Fields 




g <T d,A d v A -~V=g <T d,C d v C 



gRC 2 - l -KH 2 ^g~A 2 C 




(153) 
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The equations of motion are: 
1 5S 



A>» - kH 2 AC = , (154) 



f=gSA 

— — - C» --RC- -kH 2 A 2 



= C 5 ^ - 2H 2 C - -kH 2 A 2 = . (155) 

As usual, the perturbative initial value solution comes from iterating the 
Yang-Feldman equations: 8 

A(t,x) = A (t,x) - kH 2 J* dt'a 3 {t') J d 3 x'G rct (t,x ;t',x' ) A{x')C{x') (156) 
C(t, x) = C {t, x) - \^H 2 f G dt ' a3 (t') J d3x ' G Zt(t, x ; t', x) A 2 (x) . (157) 



The minimally coupled free field Aq is identical to ()124|) . Its conformally 
coupled cousin is: 

Co(t,x) = {e**v(t,k) 7S + e-**v*(t,k) 7 J } , (158) 

where 7! and 7^ are canonically normalized creation and annihilation oper- 
ators and v(t, k) is the conformal mode function: 

ik 



-i exp 



v{t,k) = -f= - L f° (t)j . (159) 



^2k a(t) 

The minimally coupled retarded Green's function G ret (t, x\ t', x') is the same 
one (J49|) we have been using since Section 3. The conformally coupled re- 
tarded Green's function is: 



n / ^ H 2 { 8{H\\x-x'\\ + i - , 

v > An [ aa H\\x — x'\\ J 

• Stochastic Realization: We seek a simplification of the Yang-Feldman 
equations for stochastic fields A m (t, x) and C m (t, x) which preserves the lead- 
ing infrared logarithms. The A equation ()156|) can be reduced as in the 
previous examples: 

A m (t,x) = A m0 (t,x) - ^ /' dt' A m (t',x)C m {t',x) . (161) 



^Field arguments are sometimes compressed to 4- vector notation: (t,x) = x. 
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The free stochastic mode sum A m0 (t,x) - as in the previous examples - is 
given by (|124|) . We obtain the local stochastic field equation by differentia- 
tion: 

kH 

Leading Infrared A m (t,x) = f A (t,x) — A m (t, x) C m (t, x) 

(162) 

The stochastic source f A (t,x) is given by (|136| . 

To obtain the analogous reduction for the C equation (j!57|) note first that 
a free field does not produce infrared logarithms unless its mode function goes 
like k~2 near k = 0. The conformally coupled mode function ()159|) goes like 
k~z so its associated free field truncates to zero: 

Leading Infrared =>- Co(t,x) — > C m o(t,x)=0 . (163) 

Then, we determine the infrared limit of the conformally coupled Green's 
function by beginning with the generic expression: 

G^ ct (t,x ;t',x) = i9(t-t')x (164) 

d 3 k r 

—-- [v(t, k) v*{t\ k) - v*{t, k) v(t', k, 
Using ()159j) the long wavelength expansion of the bracketed term is: 



Ak- Ax 



v(t, k) v*(t>, k) - v*(t, k) v(t>, k) = + 0(k 2 ) . (165) 

J H \aa' z a A a ) 

Consequently, the infrared limit of the conformally coupled Green's function 
is: 

a' 3 G2»(t, x ; t', a? ) — i 6{t - t') 5 3 (x - STj - ^ . (166) 

The stochastic Yang-Feldman equation for C m could be written as: 

C m (t, x) = f dt> (- - AUt', 2) ■ (167) 

2 Jo \a a 1 J 

However, a further simplification is possible. By neglecting the slow variation 
in A 2 m [t'x) relative to the rapidly growing scale factors we can perform the 
intergation over t'\ 

Leading Infrared =^> C m (t,x) = ——A 2 R (t,x). (168) 
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Substituting (|168|) for C m in (|162|) reveals the amusing fact that A m obeys 
the same stochastic field equation as the self-interacting scalar of Section 3 
with the identification A = — ^k 2 H 2 . Because the effective 4-point coupling 
is negative this system ought to exhibit a runaway instability. 

Equations (j!62|) and (jl68|) result from applying the reduction rules I-II 
(I94H95)) of Section 3 to the full equations of motion (jl54H155|) . First apply 
rules I-II (I94H95|) to the kinetic term of the minimally coupled scalar in the 
usual way: 

A% = -A - 3HA + ^ A — -3H(A m - f A ) . (169) 

By rule I f|94|) the leading contribution to the conformally coupled kinetic 
term is undifferentiated, so by rule II (j9*3j) there is no stochastic source: 

V 2 

C^ a - 2H 2 C = -C - 3HC + -v C - 2H 2 C — ► -2H 2 C m . (170) 
^ a z 

The complete stochastic field equations are therefore: 

-3H(A m -f A ) — kH 2 A m C m = , (171) 
-2H 2 C m -~ K H 2 A m C XR = 0, (172) 

and the rules I-II (I94H95|) are again shown to give correct results. 

• Perturbative Agreement: It is simple to compute the expectation values 
of squares of the full quantum fields A and C to a few orders in perturbation 
theory. The Feynman rules and relevant diagrams are shown in Figure 2; the 
explicit form of the basic vertex is: 

Vi2, 3 (x) = -ikH 2 . (173) 

The lowest order contributions are just the coincidence limits of the two 
propagators: 

(A 2 (t,f)) vEv = iA a (x;x) + 0(k 2 ) = ^{\JV + \na} + 0(K 2 ) , (174) 
(C 2 (t,f)) vEv = iA CF (x;x) + 0(k 2 ) = + 0(k 2 ) , (175) 
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Figure 2: The Feynman rules and the diagrammatic expansion to 0(n 2 ) of the 

vacuum expectation value of the squares of the fields in the scalar model 
with constraints. Solid lines correspond to the scalar field A and segmented lines to the 
scalar field C. 
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where "UV" stands for a constant, ultraviolet divergence. These results 
are completely consistent with the perturbative solution of the stochastic 
equations (IToij) and (ITfiKJ) : 

A m {t,x) = A m0 {t,x) + ^ f dt' A* R0 {t',x) + 0(k a ) , (176) 

Iz Jo 

C 1R (t,x) = -*lAl (t,x) +0(k 3 ) . (177) 

Based on the above expansions (jl76H177j) we make the following predictions 
for the next order - that is, two loop - perturbative results: 

Km) > VEV = + + °^ > ( 178 ) 

(^f))^ = ^ln 2 a + 0(. 4 ) . (179) 

(180) 

It will be interesting to check these predictions. 



7 Epilogue 

The subject of this paper has been infrared logarithms. These are powers 
of the number of inflationary e-foldings from when the system was released 
in a prepared state. They can arise in the expectation values of operators 
in quantum field theories which contain massless, minimally coupled scalars 
or gravitons. We demonstrated that infrared logarithms derive entirely from 
long wavelength quanta. In Section 3 we were able to obtain a completely 
finite system, with a conserved stress tensor, by truncating the ultraviolet 
modes from the free field of the Yang-Feldman equation. 

Infrared logarithms are very exciting because they may signal important 
quantum effects. However, their continued growth implies the breakdown 
of perturbation theory. A natural approach for obtaining non-perturbative 
information is the leading logarithm approximation in which one attempts 
to sum the series comprised of just the leading infrared logarithms at each 
order. Starobinskh has long argued that his formalism of stochastic inflation 
[TH ITS] recovers the leading infrared logarithms for scalar fields with non- 
derivative interactions. We have confirmed this using the infrared-truncated 
Yang-Feldman equation of Section 3. The leading infrared logarithms are 
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not changed when the free field mode functions are simplified to retain only 
the first nonzero term in the long wavelength expansion and a similar simpli- 
fication is done to the retarded Green's function. The fields of the resulting 
system are commuting, even though they depend upon creation and annihi- 
lation operators, and they obey precisely Starobinskii's Langevin equation. 

Whereas there is little hope of obtaining non-perturbative solutions for 
interacting quantum field theories in 3 + 1 dimensions, classical field equa- 
tions can sometimes be solved exactly. And simple classical equations are 
especially likely to possess exact solutions. In Section 4 we reviewed the 
exact solution obtained by Starobinskii and Yokoyama ^3] for the late time 
limit of a massless, minimally coupled scalar with a quartic self-interaction. 
We also reviewed Starobinskii's technique (HUES f° r predicting the infrared 
logarithms of explicit perturbative computations, and we worked out the 
general result at orders A , A 1 and A 2 . 

The remainder of the paper concerned extending Starobinskii's formal- 
ism to more complicated theories with derivative interactions and constraints. 
We have written down a set of simple rules (|94H95|) for accomplishing this. 
By going through the same procedure of first infrared-truncating the Yang- 
Feldman equations and then making the long wavelength approximation on 
their mode functions and retarded Green's function, we showed that rules 
I-II give the correct results for a scalar model with derivative interactions 
(Section 5) and for a reasonable scalar analogue of a model with gauge con- 
straints (Section 6). We also checked our rules against explicit perturbative 
computations. 

We emphasize that just because the rules I-II (I94H95|) can be stated 
generally does not mean they apply generally. For any particular model 
the stochastic Langevin equations implied by our rules should be derived 
by infrared-truncating the Yang-Feldman equations and then imposing the 
long wavelength approximation. Specific perturbative results should also be 
checked. 

The primary motivation for developing these rules is their application to 
gravitation where some interesting perturbative results already exist El • 
We close with an amusing and thought-provoking argument which is com- 
pletely independent of these results and assumes only that a stochastic formu- 
lation of inflationary quantum gravity exists. If so, then it is straightforward 
to prove that back- reaction must become important. Suppose the null hy- 
pothesis is correct: there is no significant back-reaction at any order. In that 
case the stochastic metric field is just the inflationary background plus the 
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free field mode sum of long wavelength gravitons in the transverse-traceless 
field hJJ{t,x): 



Null Hypothesis 



9ij(t,x) 



a 2 {t) S i:j + hl T {t,x) 



(181) 



Hence, the volume element operator is: 



y/-9(t,Z) = a 3 (t) [l - \hJJ{t,x)hJJ{t,x) +0(h s ) 



(182) 



Its stochastic vacuum expectation value is simple to compute: 



>J-9(t, x) 



) 



a 3 (t) 1 



4 



GH 2 \na + 0(G 2 ) 



(183) 



VEV 



7T 



For large observation times this eventually shrinks to zero, at which point 
back-reaction must have become significant. Hence it cannot be consistent 
to ignore back-reaction! 

In addition to illustrating the power of the stochastic formalism, the ar- 
gument we have just given establishes the need to include quantum gravita- 
tional back-reaction. However, this does not necessarily mean there will be 
significant modifications to the de Sitter background. It may be instead that 
non-linear corrections to the graviton field, and to the various constrained 
fields, add up to give only a small shift that approaches a constant at late 
times. Precisely this happens to the scalar model considered in Section 4. 
However, it was obvious in this model that the classical force pushing the 
scalar down its potential would eventually compensate the tendency of in- 
flationary particle production to push it ever higher. No such mechanism 
is apparent in gravity. Indeed, the physical picture seems rather to be that 
the self-gravitation from continual production of inflationary gravitons must 
eventually bring the universe to the verge of gravitational collapse. It would 
be as premature to ignore the possibility of a significant effect from quantum 
gravity on the basis of the scalar model as it would have been to discount 
the prospects for asymptotic freedom - in a somewhat different leading log- 
arithm approximation - based on the positive beta functions of QED and c/> 4 
theory. 
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8 Appendix: Perturbative Computations 



We shall derive expression (j!46|) from Section 5. The one loop contribution 
derives from the coincidence limit (ffiJl of the A-type propagator. At two loop 
order the contribution is: 



B 2 (t,x) 



2-loop 



d D x' J-g(x') g^(x') 



d' iA A++ (x;x' 



~ iA A (x'; x') — i5Z 



x 



dLiA 



A++{x; x 



d'iA A+ 4x;x') d' u iA A+ _(x-x') \ , (184) 



where we have used the Feynman rules of Figure 1 and equations 
In ()184|) . the subscript pairs ++ and H — refer to the two possible variations 
of the interaction vertex at x' as required when computing true expectation 
values in quantum field theory (23 I2H E] • The symbol 5Z represents the 
field strength renormalization necessary to absorb the ultraviolet divergence 
coming from the coincidence limit (J29|) of the propagator at the interaction 
point x'. 

We first partially integrate both the ++ and the H — sectors of 
using the identity: 



-g(x>) g^(x') 



d' iA A (x; x') d' u iA A (x;x' 



% \\l-9(?) sTi^) ^a(x; x')& v iA A (x; x') 



- iA A (x- x') 9; y-g(x') (T{x')& v iA A (x; x' 
Because the propagators obey: 

d'A^-9(x') g^(x')d'jA A++ (x;x f ) 



d',{^-g{x') g^(x')dl tA A+ 4x; x') 
and we can trivially write: 



iS D (x-x') 
, 



(185) 

(186) 
(187) 



tA^x^d'^A^x-x') = -d'^iA^x') 



188) 
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the partial integration gives: 



B 2 (t,x) 



2-loop 



iA A (x; x) 



■— iA A (x; x) — 5Z 



+ 



(189) 



^ f d D x' y/-g(x>) gr&) & v tA A (x>; x') iA A++ (x; x') - d'^ iA A+ _(x; x') } 



Note that the only surface term is at t' — because the ++ and H — prop- 
agators cancel whenever x'v lies outside the past light-cone of x M . That 
surface term vanishes if 5Z is chosen - as it must be - to cancel the one loop 
contribution to the self- mass squared at t = 0: 



5Z 



A H n - 2 T(D-l) 

2 (i^f^xfr 



7T COt( — D 
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Then we have: 



iA A (x'; x) — 5Z 



H D ~ 2 V(D-l) , 
-A ft — - ; n, — In a 



(47T) 



r 



which vanishes at f — 0. 

The first term in (j!89|) contributes a double infrared logarithm: 



(190) 



(191) 



~ iA A (x;x) - 



Aff 4 



In 2 a + O(lna) . (192) 



To see that the second term in ()189)) also contributes a double infrared log- 
arithm first note that the only the temporal derivative survives: 



H d-2 Y{D-l) 



& v iA A {x';x') = ^jT ^^.r,:' Ha' 5°, . 



(47T) 

Hence we can partially integrate on rj': 



r 



i\H n - 2 T(D-l) 
T (4 



D 
7T) 2 



H I d D x'a! D 



-i_9 



- f {iA 



A++( x i x ) 



1 



l A 2 A 



surface term at rj' = ) + 



i\ H D - 2 T{D) 



(193) 



+- v"i 



')} 



(194) 



(47T) 



^ H 2 J d D x' a' D { iA A++ (x; x') - zA 2 A+ _(x; x')} 
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The surface term contributes no infrared logarithms. In another context we 
have previously evaluated the leading infrared logarithm contribution from 
an integral of any power of the propagator [33]. For the case involving the 
squares one finds: 



iX H D ~ 2 T(D) 



4 (4 



7T 2 



f n¥T R2 1 d ° x ' a ' D { iA2a++{x] x ' ] " iA2a+ - {x] x ' ] } 

In 2 a + OQna) . (195) 



2 

XH 4 



2 6 7T 

Consequently, the leading infrared behavior of (jl89j) in D = 4 is: 

( B 2 (t, x) = ^ [ - In 2 a + O(lna) ] , (196) 

which is precisely the term appearing in expression 
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